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ABSTRACT 

The alignment of interstellar dust grains is described by the joint distribution 
function for certain "internal" and "external" variables, where the former de- 
scribe the orientation of a grain's axes with respect to its angular momentum, J, 
and the latter describe the orientation of J relative to the interstellar magnetic 
field. I show how the large disparity between the dynamical timescales of the 
internal and external variables — which is typically 2-3 orders of magnitude — 
can be exploited to greatly simplify calculations of the required distribution. 
The method is based on an "adiabatic approximation" which closely resembles 
the Born-Oppenheimer approximation in quantum mechanics. The adiabatic ap- 
proximation prescribes an analytic distribution function for the "fast" dynamical 
variables and a simplified Fokker-Planck equation for the "slow" variables which 
can be solved straightforwardly using various techniques. These solutions are 
accurate to C(e), where e is the ratio of the fast and slow dynamical timescales. 
As a simple illustration of the method, I derive an analytic solution for the joint 
distribution established when Barnett relaxation acts in concert with gas damp- 
ing. The statistics of the analytic solution agree with the results of laborious 
numerical calculations which do not exploit the adiabatic approximation. 



Subject headings: dust — polarization — ISM: magnetic fields 



1. Introduction 

Rapid advances in polarimetry and related observations have provided tantalizing clues 
about the mysterious mechanism which aligns the interstellar dust grains (for recent reviews, 
see Aitken 1996; Hildebrand 1996; Jones 1996; Goodman 1996; and Whittet 1996; see also 
Arce et al. 1997; Chrysostomou et al. 1995; Gerakines et al. 1995; Goodman & Whittet 
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1995; Goodman et al. 1995; and Martin 1995). The renewed interest in interstellar polar- 
ization has stimulated theoretical activity on the longstanding "grain alignment problem." 
The classical alignment mechanisms of Davis & Grccnstcin (1951), Gold (1952) and Purcell 
(1979 [P79]) have been reinvestigated and modified to include effects that were unknown 
or poorly understood at the time of the original investigations (Roberge, DeGraff & Fla- 
herty 1993 [RDGF93]; Lazarian 1994 [L94], 1995a,b,c; 1997a,b; Lazarian & Roberge 1997; 
Lazarian & Draine 1997). It has been demonstrated that Gold- type alignment can occur 
if the grains rotate with superthermal kinetic energies (Lazarian 1995d) and that Gold's 
mechanism operates efficiently in nonhnear Alfven waves (Roberge & Hanany 1990; L94) 
and regions undergoing supersonic ambipolar diffusion (Roberge, Hanany & Messinger 1995; 
Desch & Roberge 1997). The possibility that anistotropic radiation fields can spin up the 
grains (Dolginov & Mytrophanov 1976; Dolginov & Silant'ev 1976) has been confirmed in 
realistic simulations of the radiative torques (Draine & Weingartner 1996) and it has been 
demonstrated (Draine & Weingartner 1997) that such torques may play a direct role in the 
alignment process. 

This is the second in a series of papers, the collective purpose of which is to provide quan- 
titative predictions on the efficiencies of different grain alignment mechanisms. 1 have shown 
elsewhere (Roberge 1996) that, although it is generally impossible to infer the alignment ef- 
ficiency directly from observations, calculations of this type can be used in conjunction with 
observations to test any alignment mechanism unambiguously. The first paper in the series 
(RDGF93) described a mathematical technique which calculates the ahgnment by solving 
the Langevin equation for Brownian rotation. However, RDGF93 assumed that a grain's 
angular momentum is aligned perfectly with its axis of largest rotational inertia (its "ma- 
jor axis of inertia") due to the rapid dissipation of rotational energy by the Barnett effect 
(P79). Unfortunately, this is not an accurate assumption if the grains rotate with thermal 
kinetic energies: In an important paper, L94 pointed out that fluctuations in the Barnett 
magnetization will prevent perfect alignment even if the Barnett timescale is many orders of 
magnitude shorter than the timescales for other processes. In principle, one could incorporate 
the Barnett fluctuations directly in the dynamical simulations. However, the large disparity 
between the Barnett timescale and the timescales associated with other processes makes such 
an approach prohibitively expensive. Here I describe an approximation which eliminates this 
problem. The technique developed in this paper will be exploited in subsequent papers in 
the series on Davis-Greenstein aligment of oblate and prolate grains (Roberge & Lazarian 
1997a,b), Gold's mechanism (Roberge & Lazarian 1997c) and superparamagnetic relaxation 
(Roberge 1997). 

The structure of this paper is as follows: In §2, the timescales for relevant dynamical 
processes are briefly reviewed and the fast and slow processes are identifled. In §3, the 
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adiabatic approximation is described and applied to develop a perturbative solution to the 
Fokker-Planck equation. The application of the perturbative approach and its accuracy arc 
illustrated with a simple example in §4. The main results are discussed in §5 and summarized 
in §6. 



2. Fast and slow variables 

Consider the polarization of radiation by an ensemble of aspherical, partially-aligned 
dust grains. Throughout this paper, I will assume for concreteness that the grains are 
oblate spheroids; the generalization to other shapes is straightforward. The alignment of an 
ensemble of spheroids is characterized completely by the "Rayleigh reduction factor," 



cos C 



2 - 1 



:2-i) 



(Greenberg 1968), in the sense that R is the only orientation-dependent quantity that appears 
in the equations of transfer for the Stokes parameters. Here C, is the angle between the grain 
symmetry axis, a, and the magnetic field, and the angle brackets denote averaging over 
the grain ensemble. 

The value of C for a particular grain can be specified by giving (i) the orientation of 
a relative to the grain angular momentum, J (see Fig. 1), and (ii) the orientation of J 
relative to B (see Fig. 2). With this choice of variables, R depends, in principle, on the joint 
distribution of four angles: 0, 9, and ip. However, the rapid processional motions of a 
about J (free precession) and of J about B (Larmor precession) insure that the distributions 
of and (f) are uniform. After writing cos^ C, in terms of /3, 0, 6, and if) and averaging the 
resulting expression over %l) and 0, one finds that 

cos^ C = ^ [ 1 - cos^ e - cos^ + 3 cos^ e cos^ /?] . (2-2) 

The goal of grain alignment theory is to predict the joint distribution of 9 and /3, and hence 
R, from a model of the grain dynamics. 

The motion of a with respect to J is determined primarily by internal relaxation due 
to the Barnett effect (P79; L94; LR97). The characteristic timescale for Barnett relaxation 
in an oblate spheroid is 

^^^■•^"^^ = VKh^{h-l)J^ ^^"^^ 

(LR97), where 77 is the magnetogyric ratio of the paramagnetic moments responsible for the 
Barnett effect, h and I±_ are the inertias for rotation with J parallel and perpendicular to 
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the symmetry axis, respectively, h = I\\/I±, V is the grain volume, and K is related to the 
imaginary part of the magnetic susceptibility at frequency u by 

K ^ (2-4) 
A typical angular momentum for thermal rotation at the gas temperature, Tg, is 



•/therm = Y -^11 ^T'g. (2-5) 

Substituting Jtherm into the equation (|2-3|) , one finds that the Barnett time is typically 



where ps is the density of the grain material, a and h are the grain semiaxes parallel and 
perpendicular to the symmetry axis, respectively, 6_5 = 6/ (lO^^cm), and similarly for the 
scaling of other variables in cgs units. 

The motion of J relative to B is determined by gas damping and other external inter- 
actions on a timescale which is typically long compared to the Barnett time. For example, 
the characteristic timescale for gas damping in a cloud of particles with mass m and number 
density n is 



where fth = \j2kT^/m and V\\ is a dimensionless coefficient that depends weakly on the grain 
shape (see Appendix A). For a gas of molecular hydrogen, 

l9 ^ ^-1 rp-1/2 



tgas = 7 X 10« p,o n-^' T^r a-5 s. (2-8) 
ingly, I will refer to 6 and as the "fast variables" and J, (3, and (p as the "slow variables." 



Comparing expressions ( p-6| ) and (p-8|) , we see that tear/^gas ~ 10 ^-10 ^ typically. Accord- 
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3. The adiabatic approximation 



3.1. The grain ahgnment problem 



Let {xi,X2,X3) denote the coordinates of J in the external frame. It will be useful to 
take {xi,X2,xs) to be Cartesian coordinates, {Jx, Jy, Jz), in some problems and spherical 
polar coordinates, (J, P, 0), in others (see Fig. 2). Let / be the joint distribution function for 
6, ip, and x [with x shorthand for {xi,X2, x^)], defined so that / dxi dx2 dx^ dO dip is the joint 
probability of finding the angular momentum in dxi dx2 dx^ and the symmetry axis oriented 
in d6 dtp. By choosing an appropriate set of dimensionless time and angular momentum units, 
it is always possible to write the Fokker-Planck equation for the stationary distribution in 
the form 



where Cmt and £ext are linear differential operators. The notation {9 \ x) means that /^jnt 
involves differentiation with respect to 9 but depends on x only as a parameter and similarly 
for £ext- Here e is the ratio of the timescale for internal relaxation to the shortest charac- 
teristic timescale for external interactions. How this all works out for a particular example 
is demonstrated in §4. 

The operator Ci^t represents the motion of 9 due to the combined effects of internal 
relaxation and external processes. [] It has the form 



where Aq {J,9) and Bgg{J,9) are the diffusion coefficients for the internal relaxation mech- 
anism (e.g., Barnett relaxation) and the unspecified terms of order e represent the external 
interactions. The diffusion coefficients for internal relaxation depend on the magnitude of J 
but not on its orientation in the external frame, as the notation suggests. The operator 



represents the motion of J due to gas-grain collisions and other external interactions. The 
coefficients and Bmn generally depend on 9 as well as x. 



^Recall that external interactions generally alter 9. For example, a gas-grain collision changes J on a 
timescale that is much shorter than the grain rotation period. The impulsive motion of J relative to the 
body axes generally implies a change in 6 (see Fig. 1). 




(3-1) 




(3-2) 




(3-3) 
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3.2. Adiabatic elimination of the fast variable 

We are interested in situations wliere e ~ 10~'^-10~^ typically (see §2). A solution of 
equation (|3-1D that is accurate to C(e) should therefore be satisfactory given the various 



other uncertainties in the problem, for example in modeling the grain shape. Now if we 
merely omit all of the terms of order e, then equation (|3- 1|) becomes 



- ^ l^o{J,9) f] + l^ [BeeiJ,9) f] = 0. (3-4) 

This is the Fokker-Planck equation for internal relaxation in an isolated grain, all of the 
external interactions having been omitted. Equation ( ^-4| ) determines the 6 dependence 
of / but contains no information about the dependence on J, which appears in eq. i W^ ) 
only as an arbitrary parameter. Of course, this is appropriate for an isolated grain, whose 
angular momentum must be specified as an initial condition. Alternatively, one may say that 
equation ( P-4| ) predicts the conditional probability distribution for 6 given J. Henceforth I 
will denote the solution of equation ( |3-4|) by /int {0 \ J) to emphasize this point. 



The solution of equation ( |3-4| ) does not depend on the particulars of the internal relax- 
ation mechanism. In a steady state, the transfer of energy by internal relaxation between 
the rotational and vibrational modes of the grain establishes an equilibrium 6 distribution 
that depends only on the laws of thermodynamics and not on the kinetics of the transfer 
mechanism (see LR97 for a detailed explanation). The equilibrium distribution is 

/int id\J) = C{J) Sin e exp [-e sin2 o) , (3-5) 

where C is a normalization constant, 

and Ts is the temperature of the solid material. Two related quantities which appear fre- 
quently in the following discussion are the conditional mean value of cos^ 6, 

/'27r rn 

x{J) = d^p de cos^ e ie\j), (3-7) 

Jo JO 



and a statistic which measures the alignment of a with respect to J, 

1 



QxiJ) = I 



xiJ) 3 



(3-8) 



In summary, expression ( |3-5| ) gives /int {0 \ J) accurate to C(e). This solution was ob- 
tained by neglecting the external interactions, and hence the motion of J, in computing the 
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6 distribution. This approach assumes, effectively, that the changes in J caused by the exter- 
nal interactions are so slow that the 6 distribution relaxes adiahatically to the solution that 
would obtain if the angular momentum coordinates were frozen at their instantaneous values. 
Henceforth I will refer to this assumption as the adiabatic approximation. Note that the 
adiabatic approximation is completely analogous to the Born-Oppenheimer approximation 
for the wavefunction of a molecule: in the Born-Oppenheimer approximation, the motions of 
the nuclei are assumed to be so slow that the electronic wavefunction relaxes adiabatically to 
the time-independent wavefunction that would obtain if the nuclear coordinates were frozen 
at their instantaneous values. 



3.3. Fokker-Planck equation for the slow variables 



Once the conditional distribution is known approximately, the "full" distribution, /, 
can be determined as follows. Introduce the function /cxt {x) defined by 



/ = /int {d I J) /cxt (a;). 



(3-9) 



with J regarded as a function of x. Clearly, the rules for computing conditional probabilities 
imply that /ext is just the distribution function for the angular momentum coordinates. The 
normalization of / and /int together imply that 



^ /ext (X) = 1, (3-10) 

consistent with this interpretation. 

The dynamical equation for /ext follows directly from equation ( |3-1| ). Substitute expres- 



sion 



into equation (p-lj ) and use the fact that /ext commutes with C\at to find 

/ext {x) Ci^t {e I x) /nt {e\J) + e Ce^t {x I e) /^t {e \ J) Ut (x) = o. (3-11) 



Now multiply equation p- 1 1\ ) by dO dip and integrate over 6 and ip- The first term on the 
left side vanishes upon integration, leaving 

r # r de £ext {X I 9) /int [9 I J) /ext {x) = 0. (3-12) 
JO JO 

and carrying out the integration, one finds 



After substituting expression ( |3-3| ) for £ext {x 
that equation ( |3-12| ) is equivalent to 



d 

dx„ 



Am (a;)/ext (x) 



+ 



2 doCfYt dojf 



Bmn i^^^ fext (^) 



(3-13) 
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where 



and 



Bran (•^) 



27r 



/ de A^{e,x)f,^,{d\j) 



2-K 



di) \ d9 Bmnie,x)^,t{0\J). 



(3-14) 



(3-15) 



Expression ( p-13|) is the main result of this paper. It is a Fokker-Planck equation for the 
angular momentum coordinates wherein the fast variable, 6, has been eliminated by the 
averaging in equations (P-14| ) and (|3-15| ). Equations (|3-14|) and ( |3-15|) say that the effective 
diffusion coefficients for the angular momentum coordinates, and Bmn, are obtained by 
averaging the ^-dependent coefficients Am and Bmn in the obvious way over the equilibrium 6 
distribution. The analysis in this section establishes that the averaging procedure is accurate 
to C(e). 



3.4. Calculating the Rayleigh reduction factor 

The calculation of the Rayleigh reduction factor simplifies in the adiabatic approxi- 
mation. After multiplying expression ( p^ ) by the joint distribution, expression ( |3-9[ ) , and 
carrying out the integration over 6 analytically, one finds that 

R = J d'xU,{x)Qx{J)Qjm, (3-16) 

where J and (3 are regarded as functions of x, Qx{J) is a known function of J (see eqs. 
13-13), and 

Q,(/3)^^(cos^/3-^). (3-17) 



4. An illustrative example 

As a simple application of the approach described in §3, consider the problem of finding 
/ext when the internal relaxation is produced by the Barnett effect and the external inter- 
actions are due entirely to gas damping. This problem was solved numerically in a study 
of Barnett relaxation by Lazarian & Roberge (1997). Here I show that all the results of 
practical interest can be obtained analytically in the adiabatic approximation. 



- 9 - 



4.1. Fokker-Planck equation for the angular momentum distribution 



Adopt a system of dimensionless variables with angular momentum measured in units of 
-'therm aud time in units of tear (<Aherm) (see eqs. p-3|| and ||2-5|| ). In this system, the diffusion 
coefficients for Barnett relaxation (see §A.l) are all of order unity and those for gas damping 
(see §A.2) are expressions of order unity times 



^Bar ("Aherm) 



(4-1) 



"gas 



(compare eq. | |A3| | with eqs. | [A14] ]-[ |ATB[ | and eq. | |A4| j with eqs. | [A18| |-| [A2I| ). This verifies 
that the dimensionless Fokker-Planck equation conforms to expression ( |3- 1| ) and hence that 
the analysis of §3 applies. 



To proceed, we wish to solve equation (|3-13| ) with the coefficients Am and Bmn appro- 
priate for gas damping. To exploit the absence of a preferred direction in the external frame, 
it is advantageous to adopt spherical polar angular momentum coordinates, (J, /9,0). After 
transforming]^ the diffusion coefficients in Appendix AQ from Cartesian to polar coordinates, 
equation ( |3-13| ) becomes^ 



1 d 

2 dJ 



BjJ /ext 



where the "polar" diffusion coefficients, 

Aj = -kJ + — 
and 

Bjj = aB^ +XB\\ 
depend on J but are independent of (3 and 0. Here 



2B\ +o [B 



B 



(4-2) 

(4-3) 
(4-4) 



k = h'-^ a + X 



(4-5) 



^The transformation of the Fokker-Planck equation from one coordinate system to another is described in 
standard texts on statistical mechanics (e.g., see Risken 1984). The procedure is straightforward but tedious 
and I will only give the relevant results here. 

^Modulo the factor e; compare eqs. (P- Ij) and (P-3|). 



have actually taken a slight liberty in writing equation (4-2). The linear coefficients Ap and are 
zero, as equation (4-2) implies, but the quadratic coefficients Bpp, B^f,^, etc are not. However, there is no 
harm in omitting the nonzero quadratic terms: they represent isotropic diffusion in the angular coordinates 
and, since the linear terms are zero, they merely insure that /ext is isotropic. There is no need to include 
the quadratic terms explicitly if we assume a priori that /oxt is isotropic. 
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and 

a = l-x (4-6) 

are functions of J, 7 depends only on the grain shape (see eq. ||A17|| ), and the dimensionless 
quantities 

and 



5± = 7(l + rs/Tg 



5|l = (1 + T,/Tg) 
depend on the grain shape and dust-to-gas temperature ratio. 



(4-7) 
(4-8) 



4.2. Closed-form solution 



Equation ( |4-2| ) is a second-order differential equation for /ext- It can be simplified by 



noting that the quantity in curly brackets is the probability current along the J direction. 
Setting the current to zero yields linear, first-order differential equation for /ext, 

Aj /ext [Bjj /ext] = 0. (4-9) 

Noting that we have defined /ext to be the probability per unit coordinate interval, dJ , we 
see that the appropriate boundary condition is 

\\m /ext = 0. (4-10) 



One can easily verify that the solution of equation (|4-9|) subject to boundary condition 
( FTO|) is 



/ext(J,/?) = CJ' ^•^^ exp[F(J)] ^ 

-DJJ 

where C is a normalization constant, the factor of sin/? is required to insure that /ext is 
isotropic, 

F{J)= I A (J') dJ\ (4-12) 



Jo is a an arbitrary constant, and 



A(J) - ^ - 7- (4-13) 



This completes the solution. The expressions for Aj and Bjj imply that /ext depends on 
just two dimensionless parameters, a/h and Tg/Tg. 
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4.3. Spherical grains {a/b = 1) 



A check on the closed-form solution is provided by the special case of spherical grains. 
For spheres, the diffusion coefficients for gas damping are independent of 6 (see §A.2), so 
Barnett relaxation can be ignored. The angular momentum distribution of spherical grains 
subject only to gas damping is Maxwellian with an effective temperature Tav = | (Tg + Tg). 
In dimensionless units, the Maxwellian distribution for spheres is 

fe^t{J,P) = C J'^ sin/3 exp 



;i + 7;/Tj 



(4-14) 



To see whether the closed-form solution reproduces expression ([4-14D , note that 7 = 
h = k = 1 and B\\ = Bj_ = (1 + Tg/Tg) for spheres. Taking Jq = in equation ( [4-12| ), we 
find that 

2 7 

A = - , (4-15) 



;i + Ts/rg) 



and 



F(j) = -jV(i + Ts/rg). 

This proves that the closed-form solution is exact for spheres. 



(4-16) 



4.4. Equal gas and dust temperatures (Tg = Tg) 

If Tg = Tg, then all of the temperatures in the system are equal and the angular mo- 
mentum distribution must be Maxwellian even if the grains are not spherical. It is possible 
to prove that the closed-form solution reduces to the Maxwellian distribution when Tg = Tg 
but the derivation is tedious. As a simpler but weaker test, one can compare the statistics of 
the closed-form solution to those of the Maxwellian. The isotropy of fcxt insures that R = 0, 
as required. Consider next the statistic 

Qx= dp dJ QxiJ)fo.Mf3). (4-17) 
JO Jo 

which measures the alignment of a with respect to J (see eq. [ P-8|| ). For a Maxwellian 
distribution, Qx = Qx,Max, where 



Q 



X,Max 



2(1 - /i 



sin-^ 



i (4-18) 



(Jones & Spitzer 1967; LR97) and h = 2/ (1 + 0? /h"^) for a homogeneous spheroid. Figure 3 
is a comparison between the values of Qx,Max given by expression ( ^-1^ ) (solid curve) and 
the Qx values predicted by the closed- form solution (open circles). 
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4.5. The general case 

The general solution is shown in Figure 4, where Qx has been plotted vs. T^/Tg for 
selected values of a/b. The solid curves are the predictions of the adiabatic approximation; 
they were calculated by evaluating the integrals in expressions ( 4-12 ) and ( 4-171) numerically, 



but otherwise they are analytic results. The symbols in Figure 4 are the results of a numerical 
calculation which did not use the abiatic approximation (LR97). It took > 90 CPU hours 
on an RS/6000 workstation to compute these numerical solutions. The symbols in Figure 4 
correspond to the case e = 0.04, that is, to a case which is close to the limit, e — > 0, where the 
adiabatic approximation becomes exact. Notice that the discrepancies between the analytic 
and numerical solutions in Figure 4 are comparable to the fluctuations in the numerical 
results. 



5. Discussion 

This paper shows that, whenever the timescale for internal relaxation is small compared 
to the timescales for external processes, one can solve the grain alignment problem without 
explicitly considering the rapid "internal" motions of the grain axes relative to J. In this 
regime, the problem reduces to solving a Fokker-Planck equation for the angular momentum 
coordinates. This approach has significant practical advantages. Simulation techniques for 
solving the Fokker-Planck equation (e.g., RDGF93) must follow the dynamics with a time 
step that is typically ~ 10^^ times the smallest timescale of the motion and average the 
motion over timescales ~ 10^ times the longest characteristic time to calulate R accurate to 
a few percent. According to §2, the characteristic timescales for the internal and external 
variables differ by a factor of ~ 10^-10^, so a straightforward approach that simulated 
the internal and external motions directly would require ~ 10^^ time steps to compute each 
R value. The adiabatic elimination of the internal variables reduces the computational 
burden by a factor of e and the solution of equation ( |3-l^j| ) by simulation techniques, though 
computationally intensive, is feasible. In the next papers in the series, we follow this approach 
to predict the efficiency of the Davis-Greenstein mechanism (Roberge & Lazarian 1997a,b), 
Gold's mechanism (Roberge & Lazarian 1997c) and superparamagnetic alignment (Roberge 
1997). 



6. Summciry 



The principal results of this paper are as follows: 

1. The adiabatic approximation greatly simplifies grain alignment calculations in regimes 
where large disparities exist between the dynamical timescales associated with dif- 
ferent processes. The adiabatic approximation prescribes an analytic solution for the 
distribution function of the fast variables. The slow variables are described by a Fokker- 
Planck equation wherein the fast variables have been eliminated by a simple averaging 
procedure. 

2. The adiabatic approximation is accurate to 0{e), where e is the ratio of the dynamical 
timescales for the fast and slow variables. 

3. The adiabatic approximation provides an efficient technique for incorporating the 
physics of internal relaxation in studies of grain alignment. In a typical application, the 
adiabatic approximation reduces the computational effort required to obtain a solution 
by a factor ~ lO^-lOl The error introduced by the approximation is typically only a 
few parts in a thousand. 

4. The approach described in this paper has been used to find an analytical solution to 
a problem solved earlier using numerical techniques (LR97). Comparing the analytic 
and numerical results illustrates the simplicity and accuracy of the method advocated 
here. It also verifies the accuracy of the numerical calculations. 



This work was partially supported by NASA grant NAGW-3001. I thank Robert Lupton 
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adiabatic elimination technique. It is a pleasure to acknowledge Dave Messinger for helpful 
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A. Diffusion coefRcients 



A.l. Barnett relaxation 



The diffusion coefficients for Barnett relaxation were derived in LR97. In ordinary 
(dimensional) units, they areQ 

(Al) 



and 



Bo 



tB^J) sin 6* 



sin 9 cos 9 

^Bar(</) 



exp(^ sin 9) / yyexp(-^ y)dy + exp{-^ cos 



(A2) 



where the hats denote dimensional quantities to avoid confusion with their dimensionless 
counterparts. In a system of units where tear (-/therm) is the unit of time and Jtherm is the 
unit of angular momentum, the dimensionless diffusion coefficients are 



A, 



—J sin cos 



and 

Be 

respectively. 



J2 



sin 9 



exp(^^sin^6') [ exp{-^'^y)dy + exp{-^ 



^ COS^ I 



(A3) 
(A4) 



A. 2. Gas damping 

The diffusion coefficients for gas damping were also derived in LR97. Relative to the 
basis of the internal frame, (jb'', xj' ^ (see Fig. 1), the mean torque has Cartesian compo- 



nents 



K = -^nmbWrJ',, (A5) 
61 \ 



Al = -^nmbWrJl (A6) 



^The careful reader may notice that expressions ( |AlD and (A2) differ from the analogous expressions in 
LR97 by factors of J and J^, respectively. This is appropriate: LR97 described Barnett relaxation in a 
coordinate system fixed to the grain, where the Barnett effect appears to move J. Appropriately, LR97 
identified Ag with the fictitious mean torque required to produce this motion. Here I describe the dynamics 
in an inertial frame, where J and 6 are independent variables. In this view, Ag is the mean rate of change 
of 9, not a torque. 
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and 



(A7) 



The dimensionless factors Fy and T±_ are functions of the eccentricity, e = y'l — aP'/b^^ with 



ril(^) = ^ { 3 + 4(1 - e^)^(e) - e"^ [l - (1 - e2)^^(e)] } 



(A8) 



and 



r±(e) = TT^ { 7 - + (1 - e^)^i7(e) + (1 - 2e^) 



32 

respectively, and 



l + e" 



l_(l_e2)2^(e)]]}, (A9) 

(AlO) 



The diffusion tensor for gas damping is diagonal in the internal frame with components 



p)b DO 

yy XXI 



and 



T 

s 

t" 

J- p- 



(All) 
(A12) 

(A13) 



To transform to dimensionless variables, multiply the mean torque by tear (■Aherm) / •Aherm 
and the diffusion tensor by ^therm/^Bar (</therm)- The dimensionless mean torque is 



and 
where 



A' = -eh^Jl 



Az — ^ Jz-i 



7^r^/r||. 

The nonzero components of the dimensionless diffusion tensor are 

SL = e7(l + T;/Tg) = eS^, 



and 



Bl = e{l + %/T,)^eB^^ 



(A14) 

(A15) 

(A16) 
(A17) 

(A18) 

(A19) 

(A20) 
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FIGURE CAPTIONS 

Fig. 1 — Angles 9 and ip specify the orientation of the symmetry axis, a, with respect to the 
angular momentum, J. The "internal frame" attached to the grain has basis , 
with parallel to a and the other basis vectors oriented as shown. 

Fig. 2 — Angles j3 and specify the orientation of J with respect to the magnetic field, 
B. The "external frame," [x,y,z), has z directed along B and the remaining basis vectors 
oriented as shown. 

Fig. 3 — Values of the statistic Qx for spherical grains plotted vs. the axis ratio. Solid 
curve: exact solution (eq. ||4-18|| ). Symbols: predictions of the adiabatic approximation. 



Fig. 4 — Values of Qx for arbitrary grain shapes and dust-to-gas temperature ratios. Results 
are shown for selected values of the axis ratio, a/h, with a/h = 0.1 (top curve, open circles), 
a/b = 0.5 (middle curve, filled circles), and a/b = 0.9 (bottom curve, triangles). Solid 
curves: analytic results obtained with the adiabatic approximation. Symbols: numerical 
results from Lazarian & Roberge (1997), obtained for the case e = 0.04 without using the 
adiabatic approximation. 




Figure 1 
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Figure 2 
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Figure 3 
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